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la the long wavelength limit the current distribution in a thin 
film superconducting strip transmission line can be described by an 
inhomogeneous Fredholm equation of the second kind. By considering a 
fluxoid conservation derivation of this equation, physical insight into 
the structure of the kernel follows naturally. An approximate analytic 
solution to the integral equation is derived for a specified range of 
geometrical parameters commonly encountered in practice. The solution 
is obtained by making use of the Liouville-Neumann method of successive 
iterations and approximating the resulting series by a series involving 
powers of a defined coupling factor. It is shown that the critical cur- 
rent of the thin film superconducting strip transmission line, based on 
the calculations in the paper and 
is underestimated by less than 5 °jo 




INTRODUCTION 


Several authors have indicated that superconductive computer component's 

which are constructed in the form of thin film strip transmission lines are ad- 
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vantageous from the standpoint of switching speed and miniaturization. 

Superconducting strip transmission lines are also useful in the transportation 

of electrical information within a superconducting computer due to inherent 
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negligible loss characteristics and high group velocity. , The latter is true 

only if the film thickness is larger than^ or comparable to^ the London pene-. 

tration depth. It has also been shown that if the film thickness is less than 

the penetration depth the group velocity is appreciably decreased, making the 
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strip line useful for de^fay line memory application. In all the above devices 
it would be useful to be able to predict the total current which can be carried 
by the strip line before it becomes normally conducting. 

Several microscopic theories have been advanced which indicate that switch- 
ing in a thin film is initiated by a critical current density. Cooper and 
Marcus have shown, independently, that the problem of thin film switching is 

complicated by the fact that the current density is not constant over the cross- 
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section of the film. Using a formal Green's function approach, and employing 
the London current-field relation, Cooper derived the general inhomogeneous 
(integral ^ Fredholm) equation of the second kind describing the current density dis- 
tribution in a single film. He also obtained a specialized equation' for the case 
in which the film thickness is less than or equal to the penetration depth - the 
thin film case. .Marcus derived an identical equation by combining the Biot-Savart 
law and London's equation, and obtained computer solutions for the thin film case. 
The above methods can be used to derive an integral equation for the case of a 
strip transmission line. Due to the complexity of the kernel, analytic solutions 


t 


to the integral equation have not been found to date for either the single film 
or strip line cases. 

It is shown below that when the integral equation is derived for the strip 
transmission line, using the concept of fluxoid conservation explicitly, certain 
useful properties of the kernel can be readily deduced. These properties allow 
the approximate evaluation of the Liouville-Neumann series for a range of geome- 
tric parameters of practical interest. The analytic solution -e-fr the current 
density distribution in the strip line exhibits the same property of current 
peaking at the film edges that Cooper and Marcus found for the single film. The 
critical current of the system can be calculated from a knowledge of the current 
density in the strip line and the critical current density obtained from micro- 
scopic theory, and/or independent experimentation. 

STATEMENT OF PROBLEM 

Consider the transmission line structure in the form of two parallel cylindri- 
cal conductors shown in Fig. 1. For the sake of generality it is assumed, at first, 
that the cross-section of each conductor is arbitrary. It is further assumed that 
the wavelength of the fields propagating along the structure is much larger than 
the transverse dimensions of the structure so that a static analysis is valid. It 
is also assumed that the current density is in the z-direction and therefore is 
only a function of x and y. The latter condition is necessary in order to satisfy 
charge conservation for this quasi-static case (v*fT ^ 0). 

The derivation of the equation describing the current distribution in the 
conductors is equally valid for the cases in which: (1) the conductors are both 
normally conducting, (2) one conductor is superconducting and the other is normally 



conducting, and (3) both conductors are superconducting. It is advantageous to 
first consider the classic case in which both cylinders are normally conducting. 

Since J * cr E in both conductors, applying Faraday’s law to the dotted con- 
tour shown in Fig. 1 yields, 

J(r)- J(o) » cvi Q ~ f H(r') • [a z x dr'} (1) 

o 

in rationalized inks units, "r is the usual position vector in the x-y plane, 

a is a unit vector in the z-direction, and d*r is a differential vector line 
z 

element of integration. H("r) is the magnetic field intensity generated by cur- 
rent elements in both conductors. For the static case, the right-hand side of 
Eq. (1) is zero and therefore the current distribution in each conductor is uni- 
form. It is interesting to note that the distribution in conductor 1 is uniform 
even if conductor 2 does not have a uniform distribution (for instance if conductor 
2 is a superconductor - this case is described below), as long as conductor 1 is 
normally conducting. 

For the case in which both conductors are perfectly conducting (cr * «>) , the 
right-hand side of Eq. (1) is still zero in the static situation. How- 
ever, at some time in the past a transient existed so that it appears that infi- 
nite current densities were generated. This unrealistic situation is alleviated 
by stipulating that during the transient and afterwards the current flows in an 
infinitely thin layer at the perfect conductor surfaces with a distribution such 
that H is time independent everywhere inside the body of the conductor. Starting 
from zero initial conditions it is evident that H is zero everywhere inside of 
the cylinders. This is guaranteed by stipulating that there is no component of 
H perpendicular to the conductor’s surface. Thus, if r gl and r g2 are the position 
vectors denoting the surface of conductors 1 and 2 ^respectively and J w (r g ) the 



current per unit width at the conductor surface, then 7 ("r ) is the solution 
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to the coupled integral equations 
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While Eqs. (2) will not be solved here, some useful information can be obtained 
from their form and is summarized below: 

(1) The magnetic flux through any imaginary surface in the interior of 
either conductor is zero. 

(2) The surface current distribution in conductor 1 is influenced by 
the cross-sectional shape of conductor 1, conductor 2, and the 
separation between them. 

(3) The surface current distribution in conductor 1 is altered if conductor 

2 does not carry the same total current as conductor 1. The distribution 
in conductor 1 is altered by the presence of conductor 2 even if con- 
ductor 2 carries zero net current. 

The static case in which cr ■ » can be simulated by a dynamic situation in which 
the skin depth is much less than the trapsverse dimensions of the conductor. 

Consider the case in which both conductors are superconducting. In this situa- 
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tion 1? « jyjr in each of the cylinders is the London penetration 

depth. Thus, instead of Eq. (1), it can be shown that 

r 

[J(r) - J(o) - p 2 f H(r') • a z x dr') ] = 0 (3) 


The quantity within the brackets in Eq. (3) is the London fluxoid. Thus Eq. (3) 



expresses the principle of fluxoid conservation. Starting from the London zero 


field initial conditions, it is seen that the fluxoid associated with this or 
any other contour within the superconducting cylinder is zero. This is the exact 
analog of the zero flux condition encountered in the case of perfect conductors. 
Conditions 2 and 3 which were stated for perfect conductors remain the same for 
superconductors except that the word "surface" must be deleted. This is apparent 
since the infinite current density situation associated with a time rate of change 
of flux no longer exists in this case. 

It is convenient at this point to introduce the vector potential defined by 
the equation 

H *= Uq 1 *7 x A (4a) 


All calculations will be carried out in the Coulomb gauge. Furthermore A is chosen 
such that it is in the z-direction and thus is only a function of x and y. Thus 
in the cylinders 
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where the integration in Eq. (4b) is over the cross-sectional areas of both con- 
ductors. Using the identity 
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and the fact that the fluxoid is zero yields 


J(r) = J(o) 
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J(r ') in |r-r'|d 2 r' (6a) 
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Eq. (6a) is identical to that derived by Cooper and Marcus. The first two terms 
on the right-hand side in Eq. (6a) are independent of r, therefore the current 
density can be expressed as 



where C is an arbitrary constant determined by the total current. 


Using the concepts introduced in the above derivation the following inter- 
esting fact is apparent: The nonuniformity of the current density in the film 

(J (r) £ J(o), r £ 0) is due to the magnetic flux crossing the plane surrounded 
by the dotted contour shown in Fig. 1. For certain cross-sectional geometries 
(such as the strip transmission line which will be discussed later) this flux 
is small if both conductors carry equal and opposite currents. Thus, in these 
cases, only small variations of the current density are expected. Furthermore, 
in the case of a strip transmission line, this concept allows certain properties 
of the kernel of the integral equation to be deduced so that a closed form solu- 
tion to Eq. (6b) can be demonstrated. 

Before proceeding to this solution, consider the case in which conductor 1 
is superconducting and conductor 2 is normally conducting. As was shown previously, 
the current density in the normal conductor is uniform, so that if 1^ ^ aclc * ^ ^ 
axe the net currents and cross-sectional areas of conductors 1 and 2 respectively, 
Eq. (6a) becomes 
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It is interesting to note that for this case the distribution in conductor 1 is 
not influenced by the presence of conductor 2 for - 0. This was not true for the 
case of two perfect conductors or two superconductors, as was shown above. 


STRIP TRANSMISSION LINE OF RECTANGULAR CROSS SECTION 


Consider the strip transmission line, shown in Fig. 2, consisting of two 
parallel, infinitely long, superconducting thin films of rectangular cross - 
section. Each is of width W and thickness d, and the two are separated by a 
distance The parameters of the system are chosen to be in a range of 
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practical interest denoted by .0001 < - < .01, (pd) < 1, and < 6. A current, 

I Q ; flows into the top film, and an equal but oppositely directed current flows 

in the bottom film. For the case of the strip transmission line system described 

above, Eq. (6b) becomes 


J (x,y) - C + Ijj: 
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J(x',y') in 


(x-x') 2 + (y-y 1 ) 2 
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Using the change of variables x* — - (x* + i) in the last integral on the right- 
hand side in Eq. (8a) , and noting that 


J(x,y) 


d ^ A 
- ~2 < k < — 


- J <-x-i,y) 


(8b) 


- 2 < x < ^ 


yields 
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J(x,y) - C + 


i. P p j(x', y ') iB i fe:0 2 .. + <J7y , ? 2 l < 

47r J m J a | (x+x'+i) + (y-y')‘ 


dx* dy' (9) 


Since Od) < 1 it will be assumed that J does not vary appreciably in the x-direction 
in the film. (This assumption will be examined more closely in the next section.) 



Therefore J(x,y) « J(o,y) = J(y) and Eq. (9) becomes, after introducing the 
dimensionless variable u = (2/W) y , 





The kernel of the integral equation is an explicit function of the thickness 
to width ratio of the films, as well as the ratio of the separation distance to 
the width. It was mentioned earlier that the nonuniformity of the current density 
is due to the magnetic flux that crosses the area bounded by the dotted contour 
shown in Fig. 1. Certainly line currents (henceforth referred to as source points 
since the analysis is two dimensional) at all points in the cross section contri- 
bute to the total flux, but due to the dimensions of the transmission line that 
are being considered, contributions from anti -symmetric source points (in the two 
films) distant from the observation point tend to cancel themselves out. It 
should be expected therefore, that virtually the entire net flux will be contri- 
buted by source points that lie within some small distance from the observation 
point (x,y). This distance should be on the order of t y the separation distance 
between the two films. Thus the kernel should be sharply peaked about u » uV^and, 
based on the discussion above, it should be expected that the magnitude of the 
kernel will be a monotonically decreasing function of (|u-u f |). Since i > d and 



- *2 < x, x 1 < “ , the two dimensional kernel in Eq. (9) is negative definite. 

By considering the process by which Eq. (10) is derived from Eq. (9), it is 
apparent that the one-dimensional kernel in Eq. (10) is also negative definite. 

A plot of the magnitude of the dimensionless kernel versus (u-u*) which 
summarizes the conclusions in the above discussion, is shown in Fig. 3. This 
approximate sketch is supported by numerical analysis and shows that the kernel 
drops off to as ^ of its maximum value when lu-u 1 ] ^ i/W, and that the bulk of 
the area under the curve (95$) is contained in the range |u-u'|: 0 — 10 Note 

that approximately 50$ of the area under the curve is contained in the range 
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| u-u 1 | : 0 — — . It might be pointed out here that the approximations in this 


and other sketches in the paper will not limit the accuracy of the final analytic 
results. The sketches are given as an aid to the reader and will be useful in 
determining the errors inherent in the analytic approximations to follow. 
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Letting A = > Eq. (10a) can be rewritten as 


J(u) e C - A / J(u') 
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Equation (11) can be solved by the Liouville-Neumann method of successive itera- 
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tions: Let J v ' (u) = C. Carrying through the integration in Eq. (11) yields 

+1 i 

K(!.u-u^|) duj.}- (12a) 
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Repeating the integration in Eq. (11), but now using J' ' (u"0 yields 


J (2 >(U) =0 
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Thus, after a repetitions it can be shown that 
n +1 +1 +1 
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Therefore the current density can be expressed as 

J (u) = lim (u) 

n — oo 


(12d) 


provided the series in Eq. (12c) converges as n oo*: The range of convergence for 

this series will be discussed later. 

It is instructive to note that to zeroth order the current density, J^(u), 

is assumed to be uniform over the width of the film. The first order correction 

to this assumption, J^(u), is computed by considering the sum of the interactions 

of all the source points in the films on a particular observation point. This sum, 
+1 

A J' KCj-u-UjJ) du L , will be referred to as the coupling factor A|k(u)|. The 
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second order correction, J (u), is determined by considering the effect of all 
the source points on a particular source point before obtaining the coupling factor. 
The higher-order corrections are further expressions of the interactions of source 
points with source points. The coupling factor can be evaluated and shown to be 
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It is useful at this point to examine some of the properties of k( u) . Figure 
4 shows the absolute magnitude of the kernel plotted against u' for various values 


of u. The area under each "u" curve is |k(u)|. Since the kernel is negative 


definite, k( u) is also negative definite. Due to the narrow effective width of 
the kernel (approximately 10(— ) )> it is seen that k( u) is a weak function of 
u ^C(u) k(o^) in the central region of the film. When u is within 10 (^) of the 

edge of the film, k(u) is a strong function of u, and in fact it can be shown that 
at the edges |k(+l)| = ^ |k(o)| + e where e is a small positive number. These 
observations are summarized in Fig. 5. 

Physically this implies that since an observation point is only affected by 

those source points within a range ) > the observation points which are not 

l 

within 10 (-) of the edges are effectively in an infinitely wide film. From 

Gq. (12d) and Fig. 5 it is seen that in this region J(u) is a weak function of u 

£ 

(in an infinitely wide film J is independent of u) . However, within 10 (— ) of the 

w 

edges, the observation point is affected .by the presence of the edge of the film, 

and J(u) is a strong function of u. 

By considering a composite of Figs. 4 and 5 it is readily seen that 
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Thus in general 
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Equations (14) and (15) are very nearly true in the central portion of the film. 
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The error in the above approximations, when u is within 10 (— ) of the film edges 
will be discussed in the next section. Combining Eqs. (12d) and (15) yields 

( 16 ) 


J(u) - C jl - * K(u)| [ £ (-l) n (X|k(o)|) n l 

u n»o 



If A|*(0) |<lthe series ia A (16) converges absolutely, and the current density dis- 
tribution can be expressed in closed form as 

j<u) * c {‘ - - +$(!) i } < i7 > 

Figure 6 indicates the range for the parameters d, W, and i over which Eq. (16) 

converges. If A|k(o)|>^ the series does not converge and the Liouville-Neumann 

method is not applicable. With the help of Fig. 5 an approximate normalized 

sketch of J (u) versus u can be drawn (Fig. 7). 
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From Eq. (17) and the relation H = - p (\7 X J) > the x-component of the 
magnetic field can be shown to be 


H ( U ) - + 2 C A,[K(l-u) - K(l+u)] 


(18) 


0 W [ 1+A|k ( o) | ] 

where the plus and minus signs refer to the top and bottom films respectively. 

An approximate sketch of Eq. (18) is given in Fig. 8. 

In principle, Ampere's law can be used to evaluate H^(x,u) . However for 
the purposes of the discussion in the next section, only the form of is necessary. 
In the central portion of the films it can readily be seen that is maximum at 
x = ~ " (<& " approximately zero at x = + |, - (i + -|) . 


ERROR ANALYSIS AND CONCLUSIONS 

There are two sources of error in the analysis presented in the last section. 
The first error was introduced when it was assumed that the current density does 
not vary in the x-direction for the case ((3d) < 1 (one dimensional approximation) . 
The second error is associated with the approximations in Eqs. (14) and (15) (the 
edge approximations) . 

A. One Dimensional Approximation 

In order to evaluate this error it is convenient to use a self consistency 
argument. In other words, the one dimensional solution, Eq. (17), is resubstituted 


into the two-dimensional integral equation, Eq. (9), in order to find the x- 

direction variation of J. In practice it is more convenient to use the differ- 
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ential counterpart of Eq. (9), namely V J s P J. The x-direction variation is 
then used to find an improved y-direction variation. For the purposes of this 
paper it is not important if J is not uniform in the x-direction, as long as 
the inclusion of this variation does not appreciably affect the y-direction 
variation. It is shown below that this is, in fact, the case. 

From the discussion of the magnetic field distribution in the last section 
it is apparent that the largest x-direction variation of J occurs in the central 
portion of the film. This then is the worst case region as far as an error in 
the one-dimensional approximation is concerned. Using the self-consistent pro- 
cedure outlined above, it is apparent that in the central region of the top 
film (Hy(x ■ + ^) « 0) 

. cosh f3(x - d/2) 

J(x,y) « J(x) * J(o) (19) 

cosh p(d/2) 


d d 

It is seen that for ({3d) = 1, J(-^) « (1.5) J(+2). Thus there can be an appreci 
able x-direction variation. However, when Eq. (19) is substituted into Eq. (9) 
the improved y-variation integral equation is 
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where it can be shown that , 

i < a < (20b) 

' P(d/2) 

Thus, at worst, a » 1.04 (for the case £d » 1) in the central region of the film. 
If A — ah the magnitude of J(u) in the central region is less than 2$, smaller 
than that predicted by Eq. (17). As was pointed out previously, this error is 
less near the film edges. 


B. Edge Approximation 

In order to evaluate the error in Eqs. (14) and (15) it is convenient to 

note that due to the properties of the kernel^ which were discussed earlier^* 3 ' 

+1 ' 

K(u’) K( | u-u' | ) du 1 = K (o) J K (u) [1 - y(u)] (21) 
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where y(u)> 0. By considering more carefully the composite of Figs. 4 and 5 used 

to derive Eq. (14), it is seen that 7(u) is largest when u = ± 1, in which case 
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7(±1) =» .25. In the central portion of the film 7(u) «0. Thus it is seen that 
at the film edges Eq. (14) is not very accurate. However, if Eqs. (12d) and (21) 
are used to derive a new expression for J(±l) it is found that 

J(±.l) « C jl-A |*(1)| 1 - (.75) A |k(o) | + (.69) A 2 |*(o)| 2 - j j (22) 

Thus J(±l) is less than 5$ smaller than that predicted by Eq. (17) for all param- 
eter values under consideration. This error is less near the center of the film. 

If it is assumed that the superconducting films switch to the normally con- 
ducting state when J at any point in the films exceeds a critical value J c , it is 
apparent that the switching is initiated at the film edges. From the error 
analysis it is clear that if the J(u) given in Eq. (17) is used to calculate a 
critical current I (in terms of J ) , this value for the critical current will 
be smaller than the true critical current by less than 5$. It is noted that in 
calculating I £ the x-direction variation of J in the central portion of the film 
should be taken into account. This is easily done by combining Eqs. (17) and (19) 
as is done in the sample critical current calculation shown in the Appendix. 
can be determined by making use of a microscopic theory, 2 and/or by experimentally 
determining I for a particular choice of parameters which lie in the range defined 

above. Once J is known, I can be determined, using the above calculations for 
c c 

any other set of parameters which lie in the range under consideration in this 



paper. 
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APPENDIX 


The critical current, I c > of a superconducting strip transmission line, for 
the case d/w « 1, can easily be calculated. For this case, the current is uni- 
form over a very large percentage of the film so that from Eq. (19) it can be 
shown that 

I « J(o) WO' 1 ) 2 sinh Od/2) (A-l) 

The dimensionless variable u will not be used in this section. From Eq. (17), 

| K C'-t v ^b,)J 

and making use of the fact that \ x ffi - /2)~ [=s ^ |k(o)|, it can be seen that 

J (y = ± f ) i 

J( -- Z - « 1 + ^ A \k(o)\ (A-2) 

where the variation of J in the x-direction at the film edges has been neglected 
in line with the discussion in the last section. According to the critical cur- 
rent density hypothesis I - X when J (± W/ 2) « J . Using this hypothesis and 

o c c 

combining Eqs. (A-l) and (A-2) yields 

I c - J c WO" 1 ) 2 sinh Od/2^/l + ~ A |*C(o)| (A-3) 




The numerator of the right-hand side of Eq. (A-3) is the critical current, I , 
calculated on the assumption that the current density is uniform in the y-direction. 
It can be seen from Eq. (13) that 


K(O) 


SML 2 f A . 1 

J d/w « 1 2 


(A-4) 


From Fig. 6 it is evident that Eq. (A-4)* is satisfied for the case d/w < .001. 
It is interesting to note that in this range Ak( o) is independent of (d/w). 
Therefore, 


I c 1 . Qd) 2 z 

I : 4 v d 

cu i ^ 


-1 


(A-5) 



For situations in which (d/W) is not much less than unity, Eqs. (A-l) and (A-4) 
are not valid. In these cases the general expressions describing the variation 
of current density in the films can be used to calculate the critical current. 
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